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We investigate the locality or nonlocality of the energy transfer and of the spectral interactions 
involved in the cascade for decaying magnetohydrodynamic (MHD) flows in the presence of a uniform 
magnetic field B at various intensities. The results are based on a detailed analysis of three- 
dimensional numerical fiows at moderate Reynold numbers. The energy transfer functions, as well 
as the global and partial fluxes, are examined by means of different geometrical wavenumber shells. 
On the one hand, the transfer functions of the two conserved Elsasser energies E+ and E' are found 
local in both the directions parallel (fcy -direction) and perpendicular (fcx-direction) to the magnetic 
guide-field, whatever the B-strength. On the other hand, from the flux analysis, the interactions 
between the two counterpropagating Elsasser waves become nonlocal. Indeed, as the B-intensity 
is increased, local interactions are strongly decreased and the interactions with small fcy modes 
dominate the cascade. Most of the energy flux in the fcx-direction is due to modes in the plane at 
fc|| — 0, while the weaker cascade in the fc||-direction is due to the modes with fc|| — 1. The stronger 
magnetized flows tends thus to get closer to the weak turbulence limit where the three- wave resonant 
interactions are dominating. Hence, the transition from the strong to the weak turbulence regime 
occurs by reducing the number of effective modes in the energy cascade. 

PACS numbers: 47.27.ck, 47.65.-d, 47.35. Tv 



I. INTRODUCTION 

The existence of magnetic fields is known in many as- 
trophysical objects such as interstellar medium, galaxies, 
accretion discs, star and planet interiors, or solar wind 
(see e.g. In most of these systems, the magnetic 

fields are strong enough to play a significant dynamical 
role. The involved kinetic and magnetic Reynolds num- 
bers in these astrophysical bodies are large enough so 
that the fiows exhibit a turbulent behavior with a large 
continuous range of excited scales, from the largest where 
energy is injected towards the finest where energy is dis- 
sipated. In many cases, a strong large-scale magnetic 
field is present and induces dynamic anisotropy. Direct 
numerical simulations that examine in detail the turbu- 
lent processes in geo- and astrophysical plasmas are very 
difficult to achieve, for only rather modest Reynolds num- 
bers can be reached with nowadays computers. One way 
around this difficulty is to model the small spatial and 
temporal scales to reproduce the large-scale behavior of 
turbulent flows. A more basic understanding of turbu- 
lence is thus needed to adequately model the fiows, in 
particular when a uniform magnetic field, constant both 
in space and time, is applied. 

As a first approximation, the incompressible magneto- 
hydrodynamic (MHD) equations can be used to describe 
the evolution of both velocity, u, and magnetic field, b, 
fluctuations. In the presence of an uniform magnetic fleld 
B (magnetic fields are here expressed in velocity units), 
the Elsasser formulation of the MHD equations, with con- 
stant unit mass density, reads 

5tz± = ±B • Vz=^ - z=F . Vz± - VP + i^V^z^^ (1) 

together with V • z^ =0, where z^ = u ± b are the 
Elsasser fluctuations, and P is the total (kinetic plus 
magnetic) pressure. We assume here equal molecular vis- 



cosity v and magnetic diffusivity rj, in other words a unit 
magnetic Prandtl number (P,. = v/rj = 1). Hereafter, 
the direction along the B magnetic field is referred to as 
the parallel direction and the projection of the wavevec- 
tors along this direction is denoted fc|| , while the two di- 
rections of planes perpendicular to B are referred to as 
perpendicular directions, the wavevector projection onto 
such planes being denoted kj_ with norm fc_L = |k^|. 

For periodic boundary conditions, equations ([T|) have 
two independent invariants in the absence of molecular 
viscosity and magnetic diffusivity, namely the Elsasser 
energies : 

^^ = \\ ^^'W '^^'^ (2) 

However, when very small viscosity and magnetic diffu- 
sivity are present, it is expected that the nonlinear terms 
cascade the energies between scales, in the so-called iner- 
tial range, up to smallest ones where dissipation becomes 
effective and removes energy from the system. The rate 
at which large scales lose energy is then controlled by 
the nonlinear terms z^ • Vz^ — VP that are responsi- 
ble for coupling different scales and cascading the energy 
towards smaller and smaller scales. The nature of the 
interactions among various scales in turbulent flows that 
lead to this cascade is a long standing problem. Under- 
standing the involved mechanisms is very important to 
predict evolution of the large-scale flow behavior, and to 
estimate global quantities in astrophysical systems, such 
as the transport of angular momentum, and accretion 
rates in accretion discs. 

High-Reynolds-number hydrodynamic turbulence, of- 
ten investigated in the framework of statistically homo- 
geneous and isotropic turbulence (which can be question- 
able in natural flows), is described to first order by the 
Kolmogorov theory [3| . In this phenomenological descrip- 
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tion, interactions between eddies of similar size give the 
dominant contribution to the energy cascade. This as- 
sumption leads to an energy spectrum in k~^/'^ and an 
energy cascade rate proportional to where Urms 

is the root-mean-square of the velocity at large scale and 
L is the typical (large) flow scale. 

The cascade in MHD turbulence is more complex, es- 
pecially in the presence of a background magnetic field. 
Even in the simplest case of zero or small intensity of the 
B-field, so that isotropy could be recovered, whether the 
MHD energy cascade can be described by a phenomenol- 
ogy d la Kolmogorov is still an open question. In particu- 
lar, the assumption that interactions between similar size 
eddies (local interactions) are responsible for the cascade 
of energy to smaller scales has been questioned in turbu- 
lent MHD flows both by theoretical arguments [1, 0, Q 
and the use of numerical simulations [1, 0, Q . It has been 
shown for mechanically forced MHD turbulence there is 
a strong nonlocal coupling between the forced scales and 
the small scales of the inertial range. Moreover, the large- 
scale magnetic field generated by the dynamo action can 
also locally affect the small scales by suppressing the cas- 
cade rate in the same manner that an initially imposed 
uniform magnetic field would. In the other limit, a strong 
B-field can lead to the flow bi-dimensionnalization, with 
a drastic reduction of the nonlinear transfers along the 
uniform magnetic field. For a B-intensity (denoted B) 
well above the rms level of kinetic and magnetic fluctu- 
ations, the MHD turbulence may be dominated by the 
Alfven waves dynamics, leading to wave (or weak) tur- 
bulence where the energy transfer, stemming from three- 
wave resonant interactions, can only increase perpendic- 
ular components of the wavevectors, i.e. components in 
planes perpendicular to the B-direction (fc_L-direction), 
the nonlinear transfers along B (fc|| -direction) being com- 
pletely inhibited [§, How MHD turbulence moves 
from the weak turbulence limit, B ^ Urms, to the strong 
turbulence limit, B ^ u^ms and _B ~ (where isotropy 
could be recovered), is an open question. 

Various authors have tried to give a physical descrip- 
tion of the strong turbulence regime with B Urms- 
Iroshnikov [llj and Kraichnan [l^ first proposed a phe- 
nomenological description that takes into account the ef- 
fect of a large-scale magnetic field by reducing the rate of 
the cascade due to the short time duration of individual 
collisions of wave packets. The resulting 1-D energy 
spectrum is then given by E{k) ^ (eBY^^k^^^^ . How- 
ever, this description assumes isotropy and, while the ef- 
fect of the large-scale field is taken into account by reduc- 
ing the effective amplitude of the interactions, the inter- 
actions themselves are considered to be local. In order to 
take into account anisotropy in strong turbulence, a scale 
dependent anisotropy has been proposed (T3| . the tur- 
bulent Z;^-eddies being such that the associated Alfven 
~ l\\/B and nonlinear l±/z times are equal (the 

so called critical balance), where Z|| and l± are the typi- 
cal length scales respectively parallel and perpendicular 
to the mean magnetic field. Repeating the Kolmogorov 



arguments, one ends up with a E{k\\,kj_) ~ fcj^ en- 
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ergy spectrum with fc|| ^ fcj^ . Recently, this result has 
been generalized in an attempt to model MHD turbu- 
lence both in the weak and the strong limit, the ratio of 
the two time scales t. /r „ , being kept fixed but not nec- 
essarily of order one |14| . In an other approach to obtain 
the transition from the strong to the weak turbulence 
limit suggested time scale for the energy cascade 

is given by the inverse average between the Alfven and 



the nonlinear time scale r 



AU these 



models however assume locality of interactions that are 
also in question in anisotropic MHD turbulence A 
nonlocal model for anisotropic turbulence has been re- 
cently proposed by one of the authors [l^; it assumes 
that the energy cascade is due to interactions between 
eddies with different parallel sizes and similar perpendic- 
ular scales, while a non-universal behavior is expected for 
moderate Reynold numbers. 

Although very useful in getting a first order under- 
standing of the processes involved in a turbulent cascade, 
cascade-energy models have to be unavoidably based on 
assumptions that need to be tested. To this respect, nu- 
merical simulations of the MHD equations are very valu- 
able because they provide information about the evolu- 
tion of the fields in the whole space, something not easily 
obtained from observations. Many numerical investiga- 
tions have been performed du ring the last two decades 
[H, HO, [m, m, m, m, m, hi and, at the achieved 
Reynolds numbers, they have demonstrated that differ- 
ent power-law exponents are obtained depending on the 
chosen forcing. In this work, we use the results of numer- 
ical simulations of free decaying MHD flows at moderate 
Reynolds number to investigate the MHD interactions 
for various intensities of the external magnetic field. In 
particular, we try to investigate whether the transfer of 
energy in the parallel and perpendicular direction is lo- 
cal (i.e. the two energies cascade between nearby 
wavenumbers) or nonlocal (i.e. distant wavcnumbers are 
involved in the cascade), and whether the coupling be- 
tween the two oppositely moving waves z"*" and z~ (that 
do not exchange energy) is local or not; and if not, which 
modes are responsible for the energy cascade. 

The paper is organized as follows. In the next section, 
we give the precise definitions of the transfer functions 
and partial fluxes used to analyze the nature of the en- 
ergy cascade. The details of the numerical simulations 
are given in Section Ilia. In Section Illb, we investigate 
the locality or nonlocality of the energy transfers, and in 
section IIIc, we examine the nature of the interactions 
between the two z+ and z~ fields. We conclude and dis- 
cuss our results in Section IV. 



II. DEFINITIONS 

Our goal is to investigate the interactions among dif- 
ferent scales. To define the notion of "scale" , we use the 
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field Fourier transforms : 

^^(k)-^|z±(x)e-'^-dx3 (3) 

defined in a 27r-periodic cube, sucli as 

z±(x)=^2±(k)e"'-, (4) 

k 

Similar size eddies will be considered as the ones whose 
Fourier transform contains similar wavenumbers. 

In any basic flow interaction, three wavevectors are in- 
volved. For example, the evolution of a given Fourier am- 
plitude z+ (k) will be coupled to a z~ (p) one and cascade 
the energy to the mode z+ (q) such that the wavevectors 
satisfy k -I- p -I- q = 0. Note that the mode z~(p) does 
not gain or lose energy from this interaction since the 
two energies i?^ and E~ are separately conserved. To 
obtain the cascade mean rate, one needs to average over 
all possible triadic interactions. To get a phcnomcno- 
logical understanding of the processes at play in MHD 
turbulence, we need to know if : i) most of the ener- 
getic exchanges occur between wavenumbers such that 
|k| ~ |q| and ii) the energy flux is a result of spectral in- 
teractions of the two fields with similar wavenumbers 
or not (|k| ~ |p|). 

To address these questions, let us consider a partition 
of the wavevectors into non-overlapping sets such that 
= Uk=i'S'k = For example 5± could be the 
spherical shells of unit width and radius K, i.e. set of 
wavevectors k that have /v < |k| < A' + l. We now define 
the filtered fields (x) so that only modes in the set 5*^ 
are kept: 

z±(x)= J2 2^(k)e^''-^ (5) 
Clearly, one gets 

z±(x)=5]z±(x). (6) 

K 

The triadic interactions among the different sets, say 
5^,5*^ and 5^, are given by: 

Tf{K, P, Q) = -J z±zj . Vz±dx3 (7) 

that express the rates at which energies are trans- 
fered from to sets due to the interactions with 
the modes belonging to set. Note that the collection 
of sets S'^ and S~ need not to be necessarily the same; 
for example, 5+ could be a collection of cylindrical shells 
while could be a collection of plane sheets. Adding 
over the index P (all sets in S^), we obtain the transfer 
functions : 

T^{K, Q) = ^ %^{K, P,Q) = - f z±zT . Vz±dx3 
p 

(8) 



that give the and E transfer rates from to 
sets due to all possible interactions. Note that the z+ 
field is not exchanging energy with the z^ field, and vice 
versa, but their interaction is responsible for the redis- 
tribution of the energy among various sets. T^{K,Q) 
can give us information about the locality or nonlocality 
of the energy transfer, i.e. whether the energy is ex- 
changed by nearby sets or long-range transfers from the 
large scales directly to the small scales are also involved. 

However, the T^{K,Q) transfer functions do not give 
us direct information on the scales of the two fields z+ 
and z~ that interact and contribute to the energy cas- 
cade. To investigate the locality on nonlocality of the 
interactions between the two Elsasser counterpropagat- 
ing waves, we introduce the partial fluxes (see [2^, |29| ) 
defined as: 

= -E^'=o/^^,-^-Vz±dx-3 (9) 

that express the flux of energy out of the outer surface 
of the shell due to the interactions with the Sp shell. 
Summation over the whole collection of sets enable 
to recover the usual definition for the global fiuxes : 

= -EK'=ol4,^^-^^^dx' (10) 

In the current work, we are going to use three differ- 
ent types of wavevector collections. We first consider 
spherical shells traditionally used in studies of isotropic 
turbulence so that a set Sj^ contains the wavevectors 
k such that K < |k| < K + 1. The second collection 
of sets are cylindrical shells along the direction of the 
guiding magnetic field. In this case, the set S"^ con- 
tains the wavevectors k such that K < k± < K + 1 
(with kj_ = y^/c^ -I- fc,^ ) . Finally, we consider planes per- 
pendicular to the B-direction, so that the set con- 
tains the wavevectors k whose fc|| -component satisfies 
K < < A' + 1 (where fc|| stands for kz). 

III. NUMERICAL RESULTS 

A. Numerical setup and initial conditions 

We integrate numerically the three-dimensional incom- 
pressible MHD equations ([T]), in a 27r-periodic box using 
a pseudo-spectral method with 256^ collocation points. 
The time marching uses an Adams-Bashforth Cranck- 
Nicholson scheme, i.e. a second-order finite-difference 
scheme. The initial kinetic and magnetic fields corre- 
spond to spectra proportional to k'^exp{—k/2)'^ for k = 
[1,8], which means a flat modal spectrum for wavevec- 
tor k up to fc = 2, to prevent any favored wavevector 
at time t = 0, and the associated kinetic and mag- 
netic energies are chosen equal, namely Ey{t = 0) = 
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FIG. 1: Temporal evolutions of the total energy (top panel) 
and of the total enstrophy (bottom panel) for the four exam- 
ined intensities of the B applied field. 



Eh( t = 0) = 1/2, as in previous numerical studies (see 
[3(| and references therein). Moreover, the correlation 
between the velocity and magnetic field fluctuations, as 
measured by the cross-correlation coefficient defined by 
2 / v(x)-b(x) dx^/{Ey+Eb), is initially less than 1%. At 
scale injection, the initial kinetic and magnetic Reynolds 
numbers are about 800 for flows aX v = -q ^ A x 10""^, 
with Urms = brms ~ 1 and an isotropic integral scale 
L ^ 2tt ] k~^Ej{k)dk/ J Ey{k)dk of about tt. The dy- 
namics of the flow is then let to freely evolve. The para- 
metric study according to the intensity of the background 
magnetic field B is performed for four different values : 
-B = 0, 1,5 and 15. All the simulations are run up to a 
computational time tmax = 15, at which the loss of the 
total energy (kinetic plus magnetic) is about 95% for the 
simulation with B ^ 0, 90% for B = 1, and 83% for the 
B = 5 and _B = 15 runs. 

Figure [T] shows the time evolution of the total en- 
ergy, E{t) = {E+{t) + E-{t))/2, and the total enstrophy, 
n{t) = 1/2 /[w2(x,t) +j^(x,t)] dx^ (where w = V x u 
stands for the vorticity field, and j = V x b for the cur- 
rent), for the four different simulations. One can note 
that the influence of the strength of the external magnetic 
field clearly slows down the flow dynamics. Note that the 
B = 5 and B = 15 flows present a very similar temporal 
behavior in energy as well as in enstrophy. The analysis 
that follows in the next subsection is based on the outputs 
of the runs at i = 4 where the spectra are fully developed 
and all runs wave roughly the same enstrophy. At this 
time, the cross-correlation coefficient, that one can also 
write {E+ -E-)/{E+ + E-), is about 3.6% and 2.5% re- 
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FIG. 2: Total energy spectra in the perpendicular (top panel) 
and the parallel (bottom) directions. 

spectively for the _B = and B = 1 flows, while it is close 
to 1.6% for B — 5 and 1% for B = 15, with thus a lesser 
increase in the stronger magnetized flows. The energy 
spectra in the perpendicular direction to the uniform B- 
fleld E{k^) = 1/2 J[z-{k)]'^ + [z+{k)]'^k^dkii and in the 
parallel one E{kn) = 1/2 J[z.-{k)]'^ + [z.+ {k)]'^k^dk± are 
shown in Figure [2] at the same time. Clearly as the mag- 
netic field intensity is increased, the spectrum in the fc||- 
direction becomes steeper. Because the planes at fcy = 
and = 1 are shown to play an important role in the 
cascade, we mention here their properties in more details. 
In absence of the applied magnetic field, the modes with 

= contain 32% of the total energy and 8% of the to- 
tal enstrophy, and the = 1 modes have 30% of the total 
energy and 10% of the total enstrophy. In the strongly 
anisotropic case, B = 15, the fcy = modes contain 55% 
of the total energy and 34% of the total enstrophy while 
the = 1 modes have 37% of the total energy and 35% 
of the total enstrophy. 

In our investigation, we focus on the cascade of the E" 
energy. The cascade has also been analyzed and it 
gives qualitatively similar results. Wc consider separately 
the cascades in the perpendicular and parallel directions 
relatively to the applied magnetic field. For this reason, 
we examine three different types of flux; (i) the flux across 
spheres of radius k = |k| that corresponds to an isotropic 
analysis, (ii) the flux across cylinders of radius k = k± 
that corresponds to the flux in the perpendicular direc- 
tion and (iii) the flux across planes located at k = \k\\\ 
that corresponds to the flux in the direction parallel to 
the B-field direction. Figure [3] shows these three fluxes, 
as a function of fc, for various B-intensities. It is clear 
that as the amplitude of the large-scale magnetic field 
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FIG. 3: Fluxes Ilg{k) across : (i) spheres (solid line), (ii) 
cylinders (dashed line), (iii) planes (doted line), for B = 
(top panel), B = 1 (mid panel) and B = 5 (bottom panel). 



is increased, the parallel flux is strongly reduced. For 
-6 = 5, this flux is reduced by more than one order of 
magnitude when compared to the case with B = 0. For 
-B = 15, the parallel flux across planes is very small and 
it even takes negative values. 



O 



m 



a 



in 
m 

o 








0.010 
0.005 
0.000 

-0.005 I 

-0.010 I 
0.010 

0.005 

0.000 

-0.005 
-0.010 I 



10 
K 



20 



FIG. 4: The transfer function 1 {K,Q) that demonstrates 
the energy exchange between cylindrical shells of radius K 
and Q. Solid lines show the diagonal K — Q. The top panel 
shows the B = case and the bottom panel the B = 5 case. 



B. Energy transfers 

We now examine the locality or nonlocality of energy 
transfers from our numerical data. For two different val- 
ues of the uniform magnetic field, namely B = and 
B = 5, Figure H] shows a shadow-graph of the trans- 
fer function T~{K,Q) between and Zq, defined in 
Eqs. (O and ([H]), for energy exchanges across cylindrical 
shells (perpendicular cascade), while Figure [5] shows the 
transfer function T~ [K, Q) for energy exchanges across 
plane sheets (parallel cascade). In all cases, the trans- 
fer is concentrated along the diagonal K — Q line. This 
indicates that the cascade happens through a local en- 
ergy exchange. Similar results are obtained from the two 
other simulations at B = 1 and i? = 15 (not shown). 
Note the highly non-linear color bar used for the parallel 
cascade in the B = 5 case. This choice is due to the 
extremely fast decrease of the amplitude of T~ {K, Q) as 
the wavenumbers K and Q become large. From Figures 
m and [51 it can be seen that most of the energy exchange 
happens close to the diagonal line {K — Q). This implies 
that waves traveling in the same direction exchange en- 
ergy between similar size wavenumbers. In the strong B 
flow, some inverse cascade is also visible in the parallel 
cascade (Fig. [5]) as indicated by the dark lines below the 
diagonal and the bright ones above the diagonal. 

To get a better understanding of the T^(/\, Q) trans- 




FIG. 5; The transfer function 1 (A', Q) that demonstrates 
energy exchanges between plane sheets located at distance K 
and Q from origin, for B = (top panel) and B — 5 (bottom 
panel) cases. Solid lines indicate the K = Q diagonal. 
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FIG. 6: The energy transfer function 1 {K,Q) for the per- 
pendicular cascade (cylinders) for Q = 10 as a function of K, 
from the runs with 5 = (sohd line), B = 1 (dashed line) and 
B = 5 (dotted line) . Amplitudes of 1 (K, Q) are normalized 
to have the same order of magnitude than in the 5 = case. 
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FIG. 7: The energy transfer function T~{K,Q) for the par- 
allel cascade (planes) for Q = 10 as a function of K, from 
the data at B = (solid line), 5=1 (dashed line) and 
5 = 5 (dotted line). Amplitudes of T~(K, Q) are normalized 
to have the same order of magnitude than in the 5 = case. 



for functions, wc look at a single wavcnumber Q. Fig- 
ure [6] displays T~ (K, Q) for the perpendicular cascade 
(cylinders) at Q = 10 as a function of K, whereas Figure 
[7] shows it for the parallel cascade (planes) at Q = 10. 
To compare the results obtained from the different B 
cases, the T~ {K,Q) amplitudes are normalized so that 
all transfers are of the same order of magnitude. Positive 
values of T~ {K, Q) imply that the shell K receives energy 
from the shell Q = 10 (Fig. [51 perpendicular case) and 
(Fig. [71 parallel case) while negative values of T~ {K, Q) 
mean that the shell K gives energy to the shell (5=10. 

For the perpendicular cascade, the shell Q = 10 re- 
ceives most energy from slightly smaller wavenumbers 
than K = 10 and it gives energy to slightly larger 
wavenumbers. This implies a locality in the energy trans- 
fer, since it is mostly the nearby cylindrical shells that 



exchange energy. The parallel cascade presents a similar 
behavior; the shell Q = 10 receives energy from slightly 
smaller wavenumbers than K ~ 10 and it gives energy to 
slightly larger wavenumbers. Note however that for the 
5 = 5 flow, there is also some trace of an inverse cascade 
(energy transfer from the wavenumber Q = 10 to the 
wavenumber K ~ 8). This local behavior has also been 
found in isotropic {B = 0) decaying MHD turbulence 
simulations [3l|. Nevertheless, we need to note that in 
forced MHD turbulence where the magnetic field is gen- 
crated by dynamo action, strong nonlocal transfers also 
exist [1] . Whether these nonlocal transfers arc present 
in the forced anisotropic regime still needs further stud- 
ies. 



C. Nonlinear interactions between z+ and z 

The analysis of the energy transfer functions has thus 
shown that the energy cascades locally. As a result, each 
Elsasser field, z"*" or z~ , exchanges energy between waves 
traveling in the same direction of similar size. Nonethe- 
less, this does not mean that interactions among op- 
positely traveling waves are local. In the limit of very 
large intensities of the background magnetic field, where 
the weak turbulence theory is valid, the energy cascade 
is due to interactions with the modes in the plane at 
= 0. Therefore, modes with > 1 interact with 
modes fcy <§: 1 to cascade the energy. To that respect, 
the interactions are nonlocal since short waves (large A:||) 
interact with long waves (small fc|| ) to cascade the energy. 
To investigate how close to the weak turbulence regime 
we are, we plot in Figure [5] the total energy flux Il~{K), 
defined in Eq. (HUj), across cylinders together with the 
partial flux 1^p^q{K); defined in Eq. ^ due to interac- 
tions with just — modes. As the strength of the 
uniform magnetic field is increased, the flux due to the 
interactions with the = modes become more and 
more dominant. In the B = 15 flow, the global and par- 
tial fluxes across cylinders become almost indistinguish- 
able suggesting that interactions with the modes in the 
plane at fcy = are responsible for the energy cascade. 
This means that the flow dynamics tends to be closer to 
a weak turbulence regime where the three- wave resonant 
interactions are dominating. 

A different behavior is obtained for the parallel en- 
ergy cascade. When a mode z~ (k) interacts with a mode 
z+(p), the z^(k) energy will move to the wavevector q 
so that the relation k -I- p -I- q = holds. If however 
p belongs to the wavevector set with p|| = 0, this rela- 
tion then reads fcy + q\\ =0 in the parallel direction, i.e. 
|fcj|| = Therefore, the energy remains in spectral 

planes located at the same distance from the origin. As 
a result, interactions with the = modes cannot con- 
tribute to the energy cascade in the parallel direction. In 
this case, the closest modes to the fcy = modes are the 
ones that gives most of the energy flux. Figure [5] shows 
the total energy fiux across planes and the partial fiux 
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FIG. 8: The total energy flux 11^ (K) (solid line) across cylin- 
ders of radius K together with the partial flux IIp^q{K) 
(dashed line) for the four different values of B, from 5 = 
(top panel) up to B = 15 (bottom panel). 
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FIG. 9: The total energy flux 11^ (K) (solid line) across planes 
at fcy = K together with the partial flux Tlp^-^{K) (dashed 
line) for the four different values of B from B = (top panel) 
up to B = 15 (bottom panel). 



only due to interactions with the modes in the plane at 
/c|| = 1 (the closest to the k^^ — plane). As the ampli- 
tude of the B-field is increased, most of the parallel flux 
comes from modes with fc|| = 1. Here, we need to note 
that the flux in the parallel direction is much noisier than 
the flux in the perpendicular direction and that it often 
presents negative values (absolute values are plotted in 
the bottom panel of Fig. [5]). A thorough analysis of the 
parallel cascade would require to average many data out- 
puts which is not possible in the case of a freely decaying 
flows. Such an analysis is left for future work. 

IV. CONCLUSION AND DISCUSSION 

In this work we examine the energy cascade and the 
interactions between different scales for freely decaying 
MHD flows in the presence of a uniform magnetic held. 
Our analysis is based on data obtained from direct nu- 
merical simulations of the MHD equations with four dif- 
ferent intensities of the applied magnetic field, in an at- 
tempt to study the transition from strong to weak turbu- 
lence limit. One clearly established result is that, as the 
strength of the uniform magnetic field is increased, the 
energy spectrum becomes anisotropic with most of the 
energy concentrated in the small fc|| wavenumbers, as al- 
ready known [27| . It is further shown that the energy flux 
in the parallel direction (relatively to the uniform mag- 
netic fleld) is also strongly suppressed when the guiding 



field in introduced. 

To investigate the locality or nonlocality of the spec- 
tral interactions, we measure the transfer functions for 
the parallel and the perpendicular cascade. The transfer 
functions in the parallel and perpendicular directions are 
found local whatever the strength of the external mag- 
netic field. As a result, the coupling between modes that 
travel in the same direction is local and the energy ex- 
change occurs between similar size eddies. This behavior 
has been shown to hold in decaying isotropic MHD tur- 
bulence simulations (with B = 0) [3l|. However, in the 
presence of a mechanical forcing, strong nonlocal inter- 
actions have been observed with a direct energy transfer 
from the forced scale to the inertial range scales @, ^ . 
If this nonlocal behavior persists in the anisotropic case 
still needs further investigations. 

The locality or nonlocality of the interactions between 
oppositely moving waves (z^ and z^), that do not ex- 
change energy, is measured by means of partial fluxes 
in the parallel and the perpendicular directions due to 
the coupling in different spectral planes. This coupling 
between oppositely propagating modes does not appear 
local. As the amplitude of the applied magnetic field is 
increased, most of the interactions occur with the fcy = 
modes that are dominant in cascading the energy. Most 
of the energy flux is thus in the perpendicular direction, 
since the = modes do not contribute to the energy 
cascade in the parallel direction. Hence, the stronger 



8 



magnetized flows tends to present a dynamics close to 
the weak turbulence limit where the three-wave resonant 
interactions are responsible for the cascade process. This 
also partly explains the similar temporal evolution in the 
B = 5 and B ~ 15 regimes (see Figure [1]) since, in both 
cases, most of the cascade is due to the fcy = modes. 

For the parallel cascade, the interactions are slightly 
different. As already said, this is due to the inability of 
the /c|| = modes to cascade the energy in the parallel 
direction, fn that case, the modes with the smallest but 
not zero fc|| (fc|| ~ 1) are the ones responsible for the 
cascade. This behavior is in qualitative agreement with 
the description of a recent phenomenological model . 
However, the lack of resolution does not allow to pursue 
a quantitative comparison. 

Finally, we would like to emphasize that we analyze 
here numerical data of freely decaying MHD flows sub- 
mitted to an external magnetic field whose amplitude 
is varied, while all the other parameters are kept un- 



changed (periodic boundary conditions, unit magnetic 
Prandtl number, initial conditions and Reynolds num- 
ber). Thus, one should be cautious in any attempt to 
generalize the obtained results, e.g. forced turbulence 
could lead to different behaviors and should be studied 
separately. The results could also be dependent on the 
kinetic Reynolds number as well on the magnetic Prandtl 
number. Furthermore the use of a refined spectral grid 
in the parallel direction, allowing the presence of more 
modes with fcy ~ 0, could alter the energy cascade. 



Acknowledgments 

This work was supported by INSU/PNST and /PCMI 
Programs and CNRS/GdR Dynamo. Computation time 
was provided by fDRfS (CNRS) Grand No. 070597, and 
SIGAMM mesocenter (OCA/University Nice-Sophia). 



[1] Ya. B. Zeldovich, A.A. Ruzmaikin and D.D. Sokoloff 
"Magnetic Fields in Astrophysics" 1990 Gordon & 
Breach Science Pub. 

A.N. Kolmogorov, Dokl. Akad. Nauk SSSR 30,299 (1941) 
T.A. Yousef, F. Rincon, and A.A. Schekochihin, J. Fluid 
Mech. 575, 111 (2007) 

M.K. Verma, Phys. Reports, 401, 229 (2004). 
M.K. Verma, A. Ayyer, A.V. Chandra, Phys. Plasmas, 
12, 82307 (2005) 

A. Alexakis, P.D. Mininni and A. Pouquet, Phys. Rev. E 
72, 046301 (2005) 

P.D. Mininni, A. Alexakis and A. Pouquet, 2005, Pliys. 
Rev. E 72, 046302 

D. Carati, O. Debliquy, B. Knaepen, B. Teaca and M.K. 
Verma J. Turb. 7, 1, (2006) 

S. Galtier, S.V. Nazarenko, A.C. Newell and A. Pouquet, 
Astrophys. J. , 564, L49 (2002) 

S. Galtier, S.V. Nazarenko, A.C. Newell and A. Pouquet, 
J. Plasma Phys., 63, 447 (2000) 
P. Iroshnikov, Soviet Astron., 7, 566 (1963) 
R. Kraichnan, Phys. Fluids, 8, 1385 (1965) 
P. Goldreich and S. Sridhar, Astrophys. J. , 438, 763 
(1995) 

S. Galtier, A. Pouquet, and A. Mangeney, Phys. Plasmas 
12, 092310 (2003) 

W.H. Matthaeus and Y. Zhou, Phys. Fluids B 1 1929 
(1989) 

Y. Zhou, W.H. Matthaeus, P. Dmitruk Rev. Mod. Phys., 
76, 1015 (2004) 



[17] A. Bhattacharjee and C.S. Ng, Astrophys. J. , 548, 318 
(2001) 



[is: 

[19 

[20 
[21 

[22; 

[23; 

[24; 

[25; 

[26; 

[27 

[28; 

[29 
[30 
[31 



A. Alexakis, larXiv:0706.08T6V l (2007) 
D. Biskamp, and W.C. Miiller, Phys. Plasmas, 7, 4889 
(2000) 

J. Cho and E.T. Vishniac, Astrophys. J. 539, 273 (2000) 
P. Dmitruk, D.O. Gomez, and W.H. Matthaeus, Phys. 
Plasmas 10, 3584 (2003) 

J. Mason, F. Cattaneo and S. Boldyrev, astro-ph 
larXiv:0706.2003V l (2007) 

J. Maron and P. Goldreich, Astrophys. J. 554, 1175 
(2001) 

W.C. Miiller, and R. Grappin, Phys. Rev. Lett. 95, 
114502 (2005) 

C.S. Ng and A. Bhattacharjee, Astrophys. J. , 465, 845 
(1996) 

C.S. Ng, A. Bhattacharjee, K. Germaschewski and S. 
Galtier, Phys. Plasmas, 10, 1954 (2003) 
S. Oughton, E.R. Priest and W.H. Mattaheus, J. Fluid 
Mech. 280, 95 (1994). 

A. Alexakis, P.D. Mininni and A. Pouquet, Phys. Rev. 
Lett. 95, 264503 (2005) 

P.D. Mininni, A. Alexakis and A. Pouquet, Phys. Rev E 
74 016303 (2006) 

L.J. Milano, W.H. Matthaeus, and P. Dmitruk, D.C. 
Montgomery Phys. Plasmas 8, 2673-2681 (2001). 
O. Debliquy, M.K. Verma and D. Carati, Phys. Plasmas 
12, 042309' (2005) 



